Optical conductivity of Bose droplets in quenched disorder at T = 0: the effects of 

structure factor. 
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The optical response of a system of two-dimensional Bose condensate droplets with an internal 
charge structure interacting with impurities is considered within the mode-coupling theory. It is 
found, that the glass transition in the system of droplets occurs at much stronger disorder as 
compared to the case of stuctureless point particles. The optical conductivity of the system is found 
to exhibit a complex oscillating structure, intimately related to the charge distribution within the 
droplets. The application of the model to the problem of the far IR optical conductivity of doped 
cuprates La2-a;Sr2,Cu04 (x < 0.26) is discussed. 

Keywords: Superfluid-glass transition, memory function, Bose droplet 



INTRODUCTION 



the dynamics of the charged stripes. 
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In a recent paper—, a simple model of the superfluid- 
glass transition in the system of 2D bosons at zero tem- 
perature have been presented. The exact analytic solu- 
tion of the self-consistent equation for the memory func- 
tion permitted to retrace the quantum phase transition 
of the Bose condensate from superfluid to the insulat- 
ing (Bose-glass) regime, as the dimensionless parameter, 
measuring the strength of the short-range random poten- 
tial of impurities (pinning centers) reaches some critical 
value. 

In the present article we extend this model to consider 
the behaviour of the disordered system of Bose droplets 
with some internal charge structure, interacting with ran- 
dom impurities. The droplet may extend over a few unit 
cells on a crystalline layer and is a region of a quasi-free 
motion of a Bose particle. It may be regarded as a po- 
tential well, where the boson is confined, and the radius 
of the droplet corresponds to the relevant localization 
length. The droplet acts as a large heavy Bose quasi- 
particle (like a large polaron). To our mind, a promising 
physical example of the droplets are the local islands of 
superconducting phase of HTSC materials above Tc (see 
e.g. Ref3.), which existence is sometimes associated with 
charged stripes^. Our approach is based on the mode- 
coupling theory of diffusion and localization^, that makes 
it possible to derive the memory function of the system, 
given the bare correlation functions and the parameters 
of the random potential. 

The main qualitative result of our analysis is that the 
internal real-space structure of the droplet strongly af- 
fects the electromagnetic response of the system and gives 
rise to the novel features in optical conductivity spectra 
as compared to the case of structureless point Bose par- 
ticles. In the next section we formulate the model and 
in Section 3 show the results. In Section 4 the applica- 
tions of the model to the experimental data on microwave 
optical response of doped cuprates are discussed. We be- 
lieve, that the collective density modes in the system of 
the Bose droplets may give rise to the features in the 
optical conductivity spectra, that are usually ascribed to 



II. THE MODEL 

We consider the two-dimensional (2D) system of heavy 
Bose particles at T = 0, randomly distributed on a back- 
ground with real static dielectric constant e. Each par- 
ticle is a droplet, carrying a charge, spread over the 
droplet with the density p(r). The interaction between 
the droplets is taken into account within the random 
phase approximation, assuming, that the system is di- 
lute enough and the droplets do not overlap. 

Our goal is to study the optical response of the sys- 
tem in presence of the random impurity potential. The 
main quantity of interest in this problem is the memory 
function M{k,z) = M' + i M", that makes it possible to 
write down the optical conductivity of the system in a 
general form: 
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where LOp is the plasma frequency, z - the complex 
(Laplace) frequency. As usual, the /c-dependence is omit- 
ted in eq.([T]), so that the memory function evaluated at 
k — Q IS oi interest. Depending on the analytic proper- 
ties of M(z), eq.dl]) yields the coherent optical response 
of a clean system, or describes the diffusive charge dy- 
namics in a dirty localized phase. The common prac- 
tical method of exploring the extended-localized phase 
transition in disordered system is provided by the mode- 
coupling theory-'-, that yields the self-consistent equa- 
tion for the memory function: 
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where 7 is a bare relaxation constant, $ is the density- 
density relaxation kernel^, U{k) is the Fourier component 
of the impurity potential, d, n, m are the dimensionality 
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of the system, the concentration and mass of the con- 
stituent charge carriers, respectively. The angular brack- 
ets (. . .) stand for the average over the configurations 
of impurities. In the case of random configurations this 
reduces merely to the factor rii , the concentration of im- 
purities. It may be seen, that in absence of impurities 
M = ij and the optical conductivity, eq.([T]), reduces to 
the simple Drude form. The approach, based on eqs. (|ll2p 
have been applied to a number of model systems, both 
classical^ and quantum Fermi^ and Bose^ fluids. 

Denoting the random impurity potential as U{r) — 
Uof(r), where f(r) is dimensionless, the properly scaled 
equation ([2]) for the memory function M(0, z) of 2D Bose 
liquid at T = K takes the form.^: 
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where v = 1 for Coulomb interparticle interaction, 
V{r) ^ 1/r, and v — 2 ior logarithmic ("2D Coulomb") 
interaction, V{r) log(r). This latter type of interaction 
is typical of the topological extended structures, such as 
the vortices in the type II superconductors. The ground- 
state properties of the 2D Bose gas with logarithmic in- 
teraction have been considered in Ref.^. All quantities 
in eq.(I3]) are dimensionless, the units of momentum and 
energy being: 
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respectively. Here to* is the effective mass of the droplet, 
Q - the charge, n2D ~ their concentration and e is the real 
static dielectric constant of the medium. The properties 
of the memory function are governed by a dimensionless 
parameter 
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where rij is the concentration of impurities. 

Taking into account the internal structure of the 
droplet, the impurity potential ip{r) should be averaged 
with its charge density p{r). Hence, the Fourier trans- 
form ip{k) is multiplied by the form-factor p(fc) of the 
droplet. Thus, new expression 



ip{k)^ ^ip{kfS{k) 



(6) 



is to be put into eq.Q instead of (fik)"^, where S{k) = 
Ipfcp is hereafter called the structure factor. The form- 
factor of an isotropic 2D droplet is calculated using the 
formula: 



Pk = 2TT I p{r)Jo{kr)rdr, 
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III. RESULTS 
A. Superfluid-glass transition 

The model under consideration allows full analytic 
treatment^ for the system of point 2D logarithmic 
droplets {S{k) = 1, v = 2) interacting with pinning cen- 
ters - the random impurities with the short-range poten- 
tial (p{k) — 1. The electromagnetic properties of such a 
model share much similarity with those of pinned charge 
density waves (CDW)2,. The basic qualitative results, 
that apply to a broad range of superfluid-glass transition 
models-'--, are as follow. 

Until A is small, the low-frequency expansion of M{z) 

is: 



AI ~ miz 
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(toi = const.), and the optical conductivity ([T]) is sin- 

gular: cr' = j^^^^ d{u!), where z = oj + Qi. Thus, the 
system is superfluid at small disorder. However, as the 
full spectral weight a'{LLj)dui is constant, the reduc- 
tion of the (5(u;)-response by a factor 1 -I- toi as compared 
to the case of free Bose-condensate entails the emergence 
of the spectral features at w > 0. The spectral weight 
of the S-peak at zero frequency drops with increasing A 
and it vanishes at a critical value A = A*. At this point, 
the memory function acquires the singularity, called the 
non-ergodicity pole, and in leading order in z: 



M : 



Too 

z 



(9) 



In this regime, the optical conductivity tends to zero 
in the limit of small z, i.e. the system is in insulat- 
ing (glassy) phase. The onset of the localization cor- 
responds to the condition: toq 0. Substituting this 
into eqs. (l3l6l9p with z tending to zero, one obtains the 
following equation for A* : 
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where J is the Bessel function of the first kind. 



The solutions of eq. pU)) for point particles [S = 1) 
are summarized in Table [D In particular, for the model 
of 2D bosons with logarithmic interaction {v — 2) in 
short-range impurity potential (lyS = we recover the 
condition A* = 1, obtained in RefjJ^. Since A* varies 
depending on the model parameters, it is natural to make 
use of a renormalized invariant constant A/ A*. 

For droplets with a non-trivial structure factor, the in- 
tegral in eq. fTUl) can be solved numerically. Whatsoever 
complicated were the structure factor S'(fc), it is anyway 
obvious, that in this case A* is always greater, than the 
corresponding value for point particles, since S{k) < 1 
for non-zero k. In other words, the droplets cannot be 
localized as easily as the point particles. If the density 
rii of impurities is high, this is an odd conclusion, since 
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FIG. 1: Left panel: The dependence of the critical disorder parameter A*, ea. (|10p . on the outer radius of the droplet for some 
inner-to-outer radius ratio a: 1 - a = 0.01, 2 - a = 0.99. Right panel: The dependence A* (a) for different rout- 1 - Tout = 2.0, 
2 - Tout = 2.5, 3 - rout = 3.0. 



TABLE I: A* for point Bose particles {S{k) — 1) derived from 
ea. (|10|l for different impurity potentials 95 and interparticle 
interactions {v — 1 - Coulomb, u = 2 - logarithmic). 
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it seems likely that a large droplet gets stuck on a rugged 
potential landscape whereas a point particle may easier 
elude the obstacles. On the other hand, given a small 
number of deep potential wells, it is physically clear, that 
the point particle is trapped first, while a large droplet is 
bound less effectively, in accord with the above conclu- 
sion. Hence, making use of the terminology adopted in 
CDW theory^, we can say, that the present model corre- 
sponds the regime of strong pinning (jii low, J7o large). 

To be specific, we consider the droplets, that have the 
form of the uniformly charged rings with inner and outer 
radii denoted ri„, ro„t, respectively. The charge density 
is: p(r) = (6'(r-ri„)-6'(ro„t-r))/S', where S" ^ ■n(r1.^^~- 
rf^) is the area of the ring, whence the form-factor ([7]) is: 

_ 2 f Ji{krout) - aJi{akrout)\ 



where a = Tin /rout- The droplet is a disk for a — 0, and 
a thin ring for a close to unity. We show at Fig.l (left 
panel) the dependence A* (rout) for some values of a as 
derived from eq. (fTO|) for Coulomb droplets (i^ = 1) in the 
d{r) impurity potential (other models, listed in Table I, 
yield very similar results). As expected, the critical level 
A* of the disorder, at which the glass transition takes 
place, increases rapidly with the radius of the ring. It is 
also instructive to consider the dependence of A* on the 
thickness a of the ring, its outer radius being fixed (Fig. 1 , 
right panel). Noteworthy is the non-trivial variation of 
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FIG. 2: The optical conductivity of logarithmic (1/ — 2) 
droplets with radius rout = 15 and different ring thickness 
parameter a for A/A* = 0.1: 1-q: = 0, 2-q = 0.9, 3 - the 
spectrum of point Bose particleai. Inset: the curve 2 together 
with the properly scaled structure factor (dashed line) in rhs. 
of eq.(0. 



A* (a) curves with rout- Two competing factors, the ring 
charge density and its spatial extension, affect the cou- 
pling with impurity. As a increases from zero. A* begins 
to grow since the charge moves on average to a greater 
distance from impurity and couples weaker with it, albeit 
the charge density increases. However, this reasoning 
may become wrong for large droplets, where the char- 
acter of A* (a) dependence changes qualitatively. The 
behaviour of the curve 3 at right panel of Fig.l means, 
that the non-ergodicity transition for thin rings (large 
a) of large enough radius may occur at lower A*, than 
for disks, i.e. these rings are more readily trapped by 
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impurities. 

B. Optical conductivity 

Now consider the effects of the structure factor ^TT\i 
on the optical conductivity spectra ([1]). Some typical ex- 
amples of <t{uj) spectrum are shown at Figs. 2, 3. These 
results are obtained numerically from eas. (|ll3l6llip given 
A/ A* < 1 (superfluid regime). The spectra display 
a broad asymmetric feature, related to the collective 
(plasmon-like) density modesii^. As compared with the 
corresponding spectra of point particles with the same 
A/A* ratio, that display a single smooth feature (dashed 
curves at Figs. 2, 3), the spectra of the droplets are slightly 
shifted to lower frequencies and exhibit a peculiar os- 
cillating structure. The oscillations become more pro- 
nounced with increasing radius of the droplets and (or) 
inner-to-outer radius ratio a. 

The origin of the oscillations can be understood as fol- 
lows. The main contribution to the integral in eq.Q 
comes from the vicinity of the wave numbers A;^ > 
(* = 1,2): 

kf + k^'" - ReC = 0, C = 2(z + M(z)), (12) 

where the integrand has the peaked structure. In partic- 
ular, as far as A <C A*, |M(z)| is small and as a crude 
zero-order approximation, we may put A/ = in rhs. of 
eq.([3]), so that Im C is infinitesimal and the peaks are 
very sharp. Then, the imaginary part of eq.Q takes the 
form: 

oo 

M" = w— jY^R,{k)S{k)S{k~h)dk, (13) 
^ ^ 

where Ri{k) denotes the regular part of the integrand 
near k = ki and the structure factor is not included in 
Ri. The spectral dependence of the memory function and 
the optical conductivity ^ are thereby governed by the 
strongly oscillating expression S{ki{Lo)). E.g., for v = 2 
eas. (|12ll3p yield 

M" = 2A^^^^5((cc>2-l)i/4) , uj>l. (14) 

It may be seen at the inset of Fig. 2, that the rhs. of 
eq. reproduces correctly the oscillations of the optical 
conductivity spectrum. 

As the amount of the disorder A/ A* increases, the os- 
cillations of optical conductivity become less pronounced. 
The point is that the imaginary part M" of the memory 
function and IvciC, increase with A so that the (5-peaks 
in eq. p^ acquire a sizeable half-width A/c. Thus, the 
contributions of the structure factor S{k) add up from 
a range of wave numbers {ki{ijj) — A /c, ki{uj) -I- A /c) and 
partially cancel each other. 

We have investigated also more complex axially sym- 
metric charge distributions. All of them reduce to a dis- 
crete or continuous series of concentric rings, considered 



above. The resulting optical spectra are found to display 
a system of sharp and flat peaks of the same origin as 
those depicted in Figs. 2, 3. 



IV. DISCUSSION 

The mode-coupling theory of the 2D weakly disor- 
dered Bose condensate have been first proposed for the 
interpretation of the optical response of doped HTSC 
oxides in Ref.^^. Despite a number of model assump- 
tions, the results obtained therein have been found to 
compare well with the properly scaled optical spectra of 
YBaaCuaOr-a. 

To our opinion, the approach considered in Refi^^ 
and in the present work deserves attention in con- 
nexion with the far-IR features in optical conductiv- 
ity spectra of underdoped HTSC's, that are presum- 
ably related to the stripe excitationsii or Wigner crys- 
tal dynamicsi2. In particular, the in-plane optical con- 
ductivity of La2_a;Sr2.Cu04 (0 < a; < 0.26) is found to 
display strong resonant structure around ~ 100 cm^^, 
that consists of a main peak and a series of less pro- 
nounced peaks at higher frequenciesii. The authors in- 
terpret these features as the signature of the optical re- 
sponse of charged stripes, extending over many (~ 10^) 
cells in the Cu-0 plane. 

These results are only a part of a great body of 
findingsi^ obtained with different techniques (neutron 
scattering^*, microwave spectroscopy^^'^^, scanning tun- 
nel microscopy^^, etc.), that give convincing evidence of 
the existence and an important role of a real-space in- 
homogeneity in the physics of doped cuprates and re- 
lated systems. The electronic inhomogeneity, considered 
mainly in the form of density waves or stripes, suggests 
the percolative nature of the transition between the insu- 
lating and superconducting regions of the Tc — x phase di- 
agram of HTSC'siiii^. The latter imphes the existence of 
the nanoscale islands of superconducting phase above Tc, 
and it was proposed, that "the superconductivity arises 
predominantly in the charged stripes (or domains)"—. 

The present model of Bose droplets may be relevant 
for the description of electromagnetic response of such a 
system of local condensate patches. The droplets emerge 
on Cu-0 planes of layered HTSC's upon doping as the 
local bosons with charge Q and effective mass to*, delo- 
calized within a region of a certain spatial extension, that 
depends on the width of the potential wells created by 
doped impurities. The droplet may be either rounded or 
"stripe" -shaped - in the latter case the droplets should 
be randomly oriented so as to restore the isotropy of the 
system. At finite frequencies the optical spectrum of the 
system consists of a broad asymmetric feature (the plas- 
mon wing) due to the collective plasmon-like excitations 
of the droplet arrays. The latter is similar to the pinned 
CDW optical response, and it was recently showr^i, that 
the optical conductivity spectrum of the 2D bosons with 
logarithmic interaction in random short-range potential 
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FIG. 3; Optical conductivity of Coulomb droplets with Tout ~ 15 and different real-space structure in superfluid regime, 
A/A* = 0.25. 1 - a = (circles), 2 - a = 0.9 (rings), 3 - point particles (no structure). Left panel: short-range impurities 
{(fi = 5{r)). Right panel: Coulomb impurities {tp = 1/r). 



(i/ — 2, ip(r) ~ S{r)) exhibits just the same frequency de- 
pendence as the CDW in strong pinning regime^. How- 
ever, as far as the Bose hquid remain confined within the 
droplets, its real-space density is essentially nonuniform, 
and the effects of a non-trivial structure factor are impor- 
tant. These manifest themselves in the strong frequency 
modulation of the plasmon wing (Figs. 2, 3), governed by 
the structure factor (eq fT4|) . 

We point out that the spectra at Fig. 3 bear a good 
deal of resemblance with microwave optical conductivity 
of lightly doped La2_a;Sra;Cu04 at x ~ 0.1, that is sup- 
posedly dominated by charged stripes excitationsii. The 
in-plane concentration of the droplets is assumed to de- 
pend on the doping concentration x. Thus, x = 0.1 corre- 
sponds to n2D = 6.25 lO^^'^ cm~^ (the area of the CUO4 
cluster being 16 A^). Taking Q = 2e, m* = 3000 Wg 
and the static dielectric constant e = 40^ (e and nie 
being the electron charge and mass, respectively), we ob- 
tain the unit of energy (eqd]) Eq — 204.25 cm~^ and the 
unit of length £ = 2'Kh/po ~ 15 A. Hence, the oscillating 
structure of the spectra at Fig. 3 falls to the energy ~ 100 
cm~ and the droplet size is ~ 225 A, or about 40 CUO4 
cells. These estimates are in agreement with the energy 
of charged stripes excitations and the mean stripe size, 
reported ii>ii. 

Assuming, that each impurity creates a droplet on 
Cu02 plane and n2D ^ it follows from eq.fO, that 
A ^ 1/x. However, it is not A itself but the ratio A/ A* 
that governs the shape of the optical conductivity spec- 



tra, where A* (eqfTOl) depends on the structure of the 
droplet and in particular, on its radius. A rigorous es- 
timate of the droplet radius could be obtained by min- 
imization of free energy and is beyond the scope of the 
present work. However, since we regard the droplet as a 
region of a coherent motion of the boson, it is plausible, 
that its radius reduces with increasing doping concen- 
tration X, when the system becomes inhomogeneous at 
smaller scale, whereby A* decreases with x. As A* varies 
strongly with the droplet radius rout (see Fig. 2), this im- 
plies, that the ratio A/A* may grow with x, although A 
and A* decrease. 

At Fig. 4 we show a series of spectra, that simulates 
the variation of the microwave optical conductivity of 
La2_a;Sra;Cu04ii' with doping. The left panel corre- 
sponds to the disks (a = 0) in short-range impurity po- 
tential, (p{k) — 1. At early stage of doping (Fig.4A) there 
exists a small number of large droplets with rout — 30.5 
and Aa/A*^ — 0.05. In these conditions the spectrum 
consists of an almost flat background and a single distinct 
feature, similarly to the optical conductivity of LSCO at 
X ^ 0.07. At Fig.4B it evolves into the wing with a 
well-pronounced structure, as discussed in previous sec- 
tion. This spectrum corresponds to x ^ 0.15 and is ob- 
tained with A^ : A'^ ~ 0.07 : 0.15 and reduced droplet 
radius rout — 15, that provides Ab/A^ — 0.2. The 
spectrum at Fig.4C corresponds to x ^ 0.26 and is ob- 
tained given Ac : Ab — 0.15 : 0.26 and still smaller 
droplets with Tout = 8, so that Ac/Aq = 0.5. It is seen. 



0.5 

0.3 

0.1 
0.5 




0.5 




Q d__i_i I I I I I I I I I I I I I 



0.5 1 1.5 2 2.5 

CO (a.u.) 



(A) 



y \ III p i ^T>-h | I I I I I I I I I I I \ \ \ [ 



1 — I I I I I I I — I I I I I I I — I I I I I I — n—r 





0.5 1 1.5 2 2.5 

CO (a.u.) 



FIG. 4: Optical conductivity of Coulomb droplets at different parameters of the disorder. Left: short-range impurities. A 
- A/A* = 0.05, B - A/A* = 0.2, C - A/A* = 0.5. Right: Coulomb impurities. A - A/A* = 0.1, B - A/A* = 0.25, C - 
A/A* = 0.5. 



that the spectrum broadens, loses its fine structure and 
its magnitude grows sizeably, in agreement with exper- 
iment. Overall, a qualitative description of all essential 
points in doping dependence of the optical conductivity 
of LSCO, reported in Refi^i, becomes possible. Similarly, 
the curves at right panel of Fig. 4 have been obtained for 
rings with a — 0.5 given the following parameters: A - 
A/ A* = 0.1 , rout = 36, B - A/ A* = 0.5, rout = 15, C 
- A/A* = 0.5, rout = 7. Again, the values of A in these 
examples meet the proportions: A"^ : A^ ~ 0.15 : 0.07 
and A^ : A^ ~ 0.26 : 0.15. 

The experimental findings of Ref.^-i have been exam- 
ined also within the diagrammatic theory of charged 
stripe&iS,. In agreement with the CDW theory by 
Fukuyama and Lee^, the authors have found the opti- 
cal conductivity of an isolated stripe to display a peak 
at a pinning frequency luq: whereas the the spectra of 
an array of interacting stripes show a second peak at a 



frequency tou, that depends on the interaction strength. 
Physically, the peak at ujq corresponds to the plasmon 
wing (that is not necessarily broad and may look like 
a single resonance, as at right panel of Fig. 3), while the 
peak at uju has no direct counterpart in our model. How- 
ever, the observed microwave spectra in cuprates gener- 
ally show a much more complex structure with a number 
of resonances. 

In conclusion, we have studied the optical response of 
the liquid of interacting 2D Bose droplets in presence of 
a random impurity potential - either Coulomb or short- 
range. Each droplet is a charged Bose particle, delocal- 
ized with a charge density p{r) in a certain spatial re- 
gion (disk or ring in the present work). As discussed in 
Sec.IIIA, the behaviour of this system is similar to that 
of CDW in strong pinning regime (a dilute distribution of 
strong pinning centers). The optical conductivity spec- 
trum of the model consists of a broad asymmetric feature 
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(plasmon wing) with a rather complex oscillating shape. 
The oscillations are intimately related to the behaviour 
of the droplet structure factor S — p{k)* p{k). The main 
qualitative conclusion to be drawn from our analysis is 
that the optical response is not only related to the pe- 
culiarities of the energy level distribution, but may also 
reflect the details of the real-space inhomogeneity of the 
system, such as the charge density within the droplets of 
a (super)conducting phase. The model can possibly de- 
scribe the microwave optical response of doped cuprates, 
where the signature of the collective mode excitations 



(density waves, stripes) are observed. 
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